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Abstract 


The ‘Schrodinger’s cat’ thought experiment highlights the counterintuitive facet of quantum theory that entanglement 
can exist between microscopic and macroscopic systems, producing a superposition of distinguishable states like the 
hctitious cat that is both alive and dead. The hallmark of entanglement is the detection of strong correlations between 
systems, for example by the violation of Bell’s inequality [1]. Using the CHSH variant [2] of the Bell test, this violation has 
been observed with photons [3, 4], atoms [5, 6], solid state spins [7], and artihcial atoms in superconducting circuits [8]. For 
larger, more distinguishable states, the conflict between quantum predictions and our classical expectations is typically 
resolved due to the rapid onset of decoherence. To investigate this reconciliation, one can employ a superposition of 
coherent states in an oscillator, known as a cat state [9]. In contrast to discrete systems, one can continuously vary 
the size of the prepared cat state and therefore its dependence on decoherence. No violation of Bell’s inequality has yet 
been observed for a system entangled with a cat state. Here we demonstrate and quantify entanglement between an 
artihcial atom and a cat state in a cavity, which we call a ‘Bell-cat’ state. We use a circuit QED [10] architecture, high- 
hdelity measurements, and real-time feedback control to violate Bell’s inequality [2] without post-selection or corrections 
for measurement inefficiencies. Furthermore, we investigate the inhuence of decoherence by continuously varying the 
size of created Bell-cat states and characterize the entangled system by joint Wigner tomography. These techniques 
provide a toolset for quantum information processing with entangled qubits and resonators [11]. While recent results have 
demonstrated a high level of control of such systems [12, 13, 14], this experiment demonstrates that information can be 
extracted efficiently and with high hdelity, a crucial requirement for quantum computing with resonators [15]. 


Quantum information processing necessitates the creation 
and detection of complex entangled states. Many implemen¬ 
tations aim to encode quantum information into a register 
of physical qubits. Alternative encoding schemes using cat 
states take advantage of a cavity resonator’s large Hilbert 
space, and allow redundant qubit encodings that simplify 
the operations needed to initialize, manipulate, and measure 
the encoded information [16, 17]. The cavity mode’s state 
can be completely described by direct measurements in the 
continuous-variable basis with the cavity state Wigner func¬ 
tion [18]. We extend this concept to express an entangled 
qubit-cavity state in what we call the joint Wigner represen¬ 
tation. We construct this representation by performing a se¬ 
quence of two subsequent quantum non-demolition (QND) 
measurements (Fig. 1), where a qubit state measurement 
is correlated with a subsequent cavity state measurement. 
When working in a cavity subspace, however, complete state 
tomography may not be required and in fact many fewer 
measurements could be used to determine a state, such as 
one with a clear mapping to single-qubit observables. By 
choosing an encoding scheme where logical states of a quan¬ 
tum bit are mapped onto a superposition of coherent states 
\P) and |—/5), we can condense the joint Wigner represen¬ 
tation down to just sixteen correlations, equivalent to a 


two-qubit measurement set. Using direct fidelity estima¬ 
tion [19, 20] and CHSH Bell test witnesses [21] within this 
logical basis, we assess the degree of entanglement of the 
state. We investigate this system’s susceptibility to deco¬ 
herence by continuously increasing the cat state amplitude 
\P\. We measure a range in which correlations surpass the 
Bell violation threshold and observe its decline due to de¬ 
coherence, benchmarking the efficiency of our encoding and 
detection schemes with cat-state qubits. 

This experiment utilizes a circuit QED architecture [10, 
22] consisting of two waveguide cavities coupled to a sin¬ 
gle transmon qubit [23, 12, 13]. One long-lived cavity (re¬ 
laxation time Ts = 55 /rs) is used for quantum informa¬ 
tion storage, while the other cavity, with fast field decay 
(relaxation time = 30 ns ), is used to realize repeated 
measurements. A transmon qubit (relaxation and decoher¬ 
ence times Ti, T 2 ~ 10 /rs) is coupled to both cavity modes 
and mediates entanglement and measurement of the storage 
cavity state. All modes have transition frequencies between 
5-8 GHz and are off-resonantly coupled (see methods for de¬ 
tails). The storage cavity and qubit mode are well described 
by the dispersive Hamiltonian: 

H/h = ujsCi^CL + {ujq — ya'^'a) |e) (e| (1) 
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Figure 1: Sequential detection for entanglement characterization, (a) A quantum circuit outlines the method 
to prepare and measure entanglement between a qubit and cavity state using sequential detection. State preparation is 
performed by first creating a product state \ 2 p) = -^{\g) + |e)) (8) \/3) with a cavity displacement of amplitude P and 

a qubit gate R\ corresponding to a f rotation around the ^-axis. A conditional gate using the dispersive interaction, 
produces the entangled state \P>-b) = “/^))* Tomography is performed by measuring an observable of both 

the qubit and cavity with sequential QND measurements. A pre-rotation Ri allows qubit detection along one of three basis 
vectors X, F, and Z. The qubit is reset and a cavity observable, the displaced photon number parity Pa^ is mapped to the 
qubit for a subsequent measurement. Sequential detections are binary results compared shot-by-shot to determine qubit- 
cavity correlations, (b) The space spanned by the superposition of quasi-orthogonal coherent states \P ), \—P) constitutes 
an encoded quantum bit in the cavity. While the cavity state can be represented by its Wigner function, this logical 
state is also described by a vector within its encoded Bloch sphere. For well-separated coherent state superpositions, the 
entangled state I^Pb) is then equivalent to a two-qubit Bell state. 


where a is the storage cavity ladder operator, |e) (e| is the 
excited state qubit projector, cjs, uOq are the storage cav¬ 
ity and qubit transition frequencies, and y is the disper¬ 
sive interaction strength between the two modes (1.4 MHz). 
This interaction creates a shift in the transition frequency 
of one mode dependent on the other’s excitation number, 
resulting in qubit-cavity entanglement [24]. As described 
in Fig. 1, the system is first prepared in a product state 
\p>) = P \ e)) ^ I/d), where \g ), |e) are the ground and 

excited states of the qubit and \P) is a coherent state of the 
cavity mode. Under the dispersive interaction we allow the 
system to evolve for a time t = creating the entangled 

I'lpn) = ^{\g,l3)+ \e,-l3)) (2) 

which we call a Bell-cat state [25, 12, 13], mirroring the form 
of a two-qubit Bell state (e.g. \p)) = -^(\gg) + |ee))). 

Correlating sequential high-fidelity measurements of the 
qubit and cavity allows state tomography of the composite 
system. We use a Josephson bifurcation amplifier [26] in a 
double-pumped configuration in combination with a disper¬ 
sive readout to perform repeated QND measurements with 
qubit detection fidelity of 98.0% at a minimum of 800 ns 
between measurements. The first measurement detects the 
qubit along one of its basis vectors {X, F, Z}. This value 
is recorded and the qubit is reset to \g) using real-time feed¬ 
back (see methods). The displaced photon-number parity 
observable of the cavity is subsequently mapped onto 
the qubit using Ramsey interferometry [18] before a second 
qubit state detection. The cavity observable Pa = DaPDj^ 
where Da is the displacement operator and P the photon 
number parity operator is detected with 95.5% fidelity (see 


methods). The Wigner function W{a) = ^ {Pa) is con¬ 
structed from an ensemble of such measurements with dif¬ 
ferent displacement amplitudes a. The correlations between 
the qubit and cavity states make up what we refer to as the 
joint Wigner functions: 

Wiia) = f (3) 

where is an observable in the qubit Pauli set 
{/, X, F, Z}. These four distributions are a complete 
representation of the combined qubit-cavity quantum state 
(see Fig. 2). While other representations exist for sim¬ 
ilar systems [27, 28, 29, 30], Wi{a) is directly measured 
with this detection scheme and does not require a den¬ 
sity matrix reconstruction. By an overlap integral (see 
methods), we determine the fidelity to a target state T = 
('0b|p|'0b) = where Wf{a) are 

the joint Wigner functions of the ideal state \Pjb) and Wi{a) 
are the measured joint Wigner functions (normalized), yield¬ 
ing a state fidelity T = (87 ± 2)% for a displacement am¬ 
plitude P = v^. This amplitude was chosen to ensure or¬ 
thogonality between logical states \ {P\—P)p = 6 x 10“^ ^ 1 
with minimal trade-off due to photon loss. Furthermore, the 
efficiency of our detection scheme can be quantified by the 
visibility of the unnormalized joint Wigner measurements 
V = ^ f (IPa}d^(y = (85 ± 1)%. Visibility V is primar¬ 
ily limited by measurement fidelity and qubit decoherence 
between detection events (see methods). The parameters 
P and V represent critical benchmarks for creating and re¬ 
trieving information from entangled states. 

The number of measurement settings required to perform 
cavity state tomography can be resource intensive. Restrict- 
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Figure 2: Joint Wigner tomography of a Bell-cat state, (a) The set of joint Wigner functions Wi{a) = - {(liPa) 
represents the state of a qubit-cavity system with correlations between the qubit = {/ , X , T , Z} and cavity 
reported for a state \'iPb) and displacement amplitude p = \/3. Shown are measurements comprised of four grids of 6500 
correlations between the qubit and cavity states. Interference fringes in (XPa) and (YPa) reveal quantum coherence in 
the entangled state, (b) A density matrix reconstruction shows the combined qubit-cavity state p in the Fock state basis, 
(c) Projecting onto the logical basis \p) {P\ -h \ —p) {—P\i this state can be further reduced exhibiting the traditional Bell 
state form. 


ing to an encoded qubit subspace, only four values of the 
cavity Wigner function W(a) are required to reconstruct the 
state, known as a direct fidelity estimation (DFE) [20, 19]. 
For large cat states | {P\—/3) p <C 1, the encoded state ob¬ 
servables map to cavity observables as: 

Xc = Po Ic = Pp p P-j5 (4) 

Yc = P jji Zc = pp — P—f3 

where {/c, Xc, Tc, ^c} form the Pauli set for the encoded 
qubit state in the cavity (see methods). Cuts in the joint 
Wigner function (Fig. 3) show these observables and their 
correlations to the qubit as a function of cat state size. As 
the superposition state is made larger, interference fringe 
oscillations increase while fringe amplitude decreases due to 
photon loss. For a state I^Pb) with \^\ = a/ 3, we estimate a 
direct fidelity JbPE = i((//c) + {XX,) - {YY,) + (ZZ,)) = 
(72 ± 2)% putting a fidelity bound on the target state with 
no corrections for visibility. This estimate is related to the 
benchmarks reported above Xdfe ~ V x X and far surpasses 
the 50% threshold for a classically correlated state. This 
indicates both high fidelity state-preparation and measure¬ 
ment, and demonstrates that strong correlations are directly 
detectable using joint Wigner tomography. 


To place a stricter bound on observed entanglement, we 
perform a CHSH Bell test on the measured state. Although 
proposed to investigate local hidden variable theory, the Bell 
test also serves to benchmark the performance of a quantum 
system that creates and measures entangled states [8, 7, 31]. 
Classical theory dictates that the sum of four correlations 
will be bounded such that: 

-2 < O = (AAe) + (AB,) - (BA,) + (BB,) < 2 (5) 

where in this experiment A, B are two qubit observables 
and Ac, Be are two cavity observables. We perform two Bell 
tests (Fig. 4) with correlations taken shot-by-shot with no 
post-selection or compensation for detector inefficiencies. In 
the first, we take observables X(^), Z(^), Xc, Zc and sweep 
both qubit detector angle 0 (see methods) and cat state 
amplitude We observe a Bell signal with a maximal vi¬ 
olation Oi = 2.30 ± 0.04 at 6> = — J for = 1. In the 
second Bell test, we follow a scheme similar to Ref. [21] and 
choose observables X, T, Xc{a)^Yc{a) where a is a displace¬ 
ment amplitude corresponding to a rotation of the encoded 
cavity state detector (see methods) and observe a maxi¬ 
mal violation O 2 = 2.14 ± 0.03 for /3 = 1. As predicted, 
a lower Bell signal is observed in the second test due to 
its greater sensitivity to photon loss, yet in both tests two 
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Figure 3: Qubit-cavity correlations, (a) Corre¬ 
lations are measured for entangled states I^Pb) with 
displacement amplitudes ranging from = 0 to 2. 
Cuts in joint Wigner functions {IPa) and {ZPa) at 
Im((a) = 0 show the increasing separation of the 
coherent state superpositions, whereas (XP^) and 
{YPa) at Re((a) = 0 reveal the interference fringe 
oscillations dependence on cat state size, (b) Sin¬ 
gle cuts dit ^ show single-shot correlations 

(crosses) with their ideal trends (solid line). Us¬ 
ing just individual measurement settings (circled), 
joint observables such as {IIc^ XXc^YYc, ZZc} of 
the qubit-cavity state can be determined. 
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Figure 4: Bell tests with a cat state, (a) Using correlations between qubit state observables X{0) and Z{0) and the 
encoded state observables Zc, Xc] a CHSH Bell test O = {AAc) + {ABc) — {BAc) + {BBc) is performed as a function of 
qubit detector angle 0. (b) Furthermore, we observe a violation with an additional Bell test using qubit observables X, 
Y and cavity state observables Xc(a), Yc{a) where a corresponds to a tomography displacement amplitude serving as a 
rotation of the effective cavity detector angle, (c-d) Both tests are performed for different cat state amplitudes P and show 
the dependence of the entangled state with photon loss and detector visibility. Squares represent measured values with 
height denoting their statistical uncertainty. Solid lines describe the predicted trends given the measured cavity decay rate 
and detection visibility. While the ideal behaviour (red) for an entangled state approaches O = 2v^, photon loss (green), 


detector visibility (blue), and their combined effects (black) 

regimes are evident. For small cat state amplitudes, the 
initial Bell signal is limited by the non-orthogonality of the 
coherent state superpositions (see methods), while for large 
displacements the system’s sensitivity to photon loss results 
in a reduction of the Bell signal. Larger, more distinguish¬ 
able states quickly devolve into a classical mixture due to 
the onset of decoherence, corresponding to the resolution of 
Schrodinger’s thought experiment. However, for intermedi¬ 
ate cat state sizes, we violate Bell’s inequality beyond the 
statistical uncertainties in both tests. 

In this letter, we have demonstrated and quantified the 
entanglement between an artificial atom and a cat state in 
a cavity mode. We determine the entangled state using se¬ 
quential detection with high-fidelity state measurement and 
real-time feedback on the quantum state. We benchmark 
the capabilities of this detection scheme with direct fidelity 
estimation and Bell test witnesses, which both reveal non- 
classical correlations of our system. This work demonstrates 
the viability of using and measuring redundantly encoded 
states in multi-level systems [16]. This implementation pro¬ 
vides a vital resource for quantum state tomography and 
quantum process tomography of continuous-variable sys¬ 
tems and creates a platform for measurement based quan¬ 
tum computation and quantum error correction using su¬ 
perconducting cavity resonators [11]. Finally, these features 


will ultimately limit the maximum Bell signal achieved. 

are directly extendible to multi-cavity systems which will 
require entanglement detection between continuous variable 
degrees of freedom and entanglement distribution of com¬ 
plex oscillator states. 
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Supplemental Materials: 

Violating Bell’s inequality with an artificial atom and a cat state in a cavity 


1 Materials and methods 

Measurement setup: Experiments are performed in a 
cryogen-free dilution refrigerator at a base temperature of 
^ 10 mK. Our output signal amplification chain consists of 
two stages. A Josephson bifurcation amplifier (JBA) [26] 
operating in a double-pumping configuration [32, 33] serves 
as the first stage, which is followed by a high electron mo¬ 
bility transistor (HEMT) amplifier. 

Eabrication techniques of the transmon qubit and the de¬ 
sign of storage and readout resonators follow the methods 
described in [ 12 ]. The refrigerator wiring (see Eig. S 2 ), in¬ 
cluding the filters and attenuators used, are similar to that 
of [13], but with the addition of a feedback system, the de¬ 
tails of which are discussed in a following section. 

Qubit-cavity parameters: The two-cavity, single-qubit 
system is well described by the approximate dispersive 
Hamiltonian: 


cavity. This allows a measurement of the time-dependent 
overlap of the cavity state with its ground state | 0 ) depen¬ 
dent on time. The lifetime of the readout cavity is found 
from its line-width. The thermal population of the qubit 
is determined from a histogram of one million single-shot 
measurements of the qubit thermal state, where the signal- 
to-noise ratio provided by the JBA allows discrimination 
between \g) and all states not \g). The thermal popula¬ 
tion of the storage cavity is found by taking the difference 
between parity measurements of the thermal and vacuum 
states of the cavity. A vacuum state is prepared by first 
performing two parity measurements on the thermal state 
and then post-selecting such that all results give even par¬ 
ity, projecting the thermal state onto |0). Einally, the known 
thermal population of the readout cavity is bounded by the 
dephasing rate T^ of the qubit: T(jy=nthK, where fith is the 
readout cavity’s thermal occupation and k is the readout 
single-photon decay rate [36]. 


Hfh = Ugolas + uOr^lar + ujqb^b 

- a/ - 

- Xqsalasb'^b - Xqrdlarb'^b - Xrsalagalar 


(SI) 


Where ujs,r,q are the storage, readout, and qubit transition 
frequencies, ag^ar^b are the associated ladder operators, and 
AT, X are the modal anharmonicities and dispersive shifts 
respectively. Table 1 details the Hamiltonian parameters 
of our system. The resonant frequency of the readout res¬ 
onator C(;^/27r is determined by transmission spectroscopy. 
The qubit frequency Ugj^'K and storage cavity frequencies 
ujql27r are found using two-tone spectroscopy. 

Qubit anharmonicity Kq is measured using two-tone spec¬ 
troscopy to observe the 0 — 2 two-photon transition [ 22 ]. 
Storage cavity anharmonicity Kg is determined by displac¬ 
ing the cavity with a coherent state and observing its time 
evolution with Wigner tomography. The resulting dynamics 
are characterized by state reconstruction and Kg is observed 
by the state’s quadratic dependence of phase on photon 
number. Einally, we predict the readout cavity anharmonic¬ 
ity Kr using its approximate dependence on the measured 
values of Kq and the qubit-readout dispersive shift Xqr [34]. 

The dispersive shift between the qubit and the readout 
resonator Xqr is found by taking the difference in frequency 
between the readout resonance when the qubit is in the 
ground and excited state. The dispersive shift between the 
qubit and the storage resonator Xqs is found using two meth¬ 
ods: photon number dependent qubit spectroscopy [35], and 
observing qubit state revival using Ramsey interferometry 
[ 12 ]. Einally, Xrs is predicted using its approximate rela¬ 
tionship between Kg and K^ [34]. 

Lifetimes and thermal populations: The lifetime of 
the storage cavity is determined by displacing to a coherent 
state, waiting a variable length of time, and then applying 
a qubit rotation conditioned on zero photons in the storage 


Measurement fidelities: We define singleshot measure¬ 
ment fidelity as Fq = ^ where P{g\g) and P{e\e) 

are the probabilities to get \g) (|e)) knowing that we start 
with 1^) (|e)). The state \g) is prepared through purifica¬ 
tion of the qubit thermal state with realtime feedback (see 
the following section). Given a preparation of |^), we have 
a 98.5% chance of measuring \g) again {P{g\g) = 0.985). 
Likewise, we find P{e\e) = 0.975 by preparing \g) and rotat¬ 
ing the state to |e). This gives a single-shot measurement 
fidelity of Fq = 98%. We find our cavity parity measure¬ 
ment fidelity by purifying the storage cavity thermal state 
into | 0 ) then performing one of two kinds of parity measure¬ 
ment (see Eig. 2 ). We report a parity measurement fidelity 
for n = 0 photons as Pc = ^^ 9 \Ei)+P{e\E 2 ) _ 95 5 %^ where 
P{g\Pi) {P{e\P 2 )) is the probability to measure \g) (|e)) 
given that the parity is even for each of the two measure¬ 
ment settings. We expect Pc to decrease with increasing 
numbers of photons in the cavity due to single photon loss 
during the measurement sequence. 

Directly from these readout fidelities, the estimated visi¬ 
bility [37] for correlated observables Vest = ( 2 Eg — l)(2Fc — 
1) = 87%. This allows us to predict the maximum Bell viola¬ 
tion possible given only measurement inefficiencies Omax = 
2 a/2 Vest = 2.47. In practice, V is directly related to the 
contrast of the joint Wigner function (see Sec. 3) which we 
measure to be 85%. This discrepancy is due to qubit de¬ 
coherence, which is studied further in Sec. 2 and puts a 
more conservative estimate for the maximum Bell violation 
achievable: Omax = 2 a/2V = 2.40. 

I/O control parameters: As shown in Eig. S3, we em¬ 
ploy a field-programmable gate array (EPGA) in order to 
implement an active feedback scheme. We use an X6-1000M 
board from Innovative Integration which contains two 1 
GS/s ADGs, two 1 GS/s DAG channels, and digital in¬ 
puts/outputs all controlled by a Xilinx VIRTEX -6 EPGA 
loaded with custom logic. We synchronize two such boards 
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in a master/slave configuration to have IQ control of both 
the qubit/storage cavity. IQ control over the readout cav¬ 
ity is performed with a Tektronix AWG, which is triggered 
by the master board. The readout and reference signals 
are routed to the ADCs on the master board, where after 
the FPGA demodulates the signal and decides whether the 
qubit is in \g) or |e). The feedback latency of the FPGA 
logic (last in, first out LIFO) is 320 ns. Additional delay for 
active feedback include cable delay (^100 ns) and readout 
pulse length with resonator decay time (320 ns). Thus, in 
total the qubit waits r^ait ^ 740 ns between the time pho¬ 
tons first enter the readout resonator and the time at which 
the feedback pulse resets the qubit. 

Implementations of feedback: Feedback is used three 
times during a single iteration of the experiment. Prior to 
the state preparation (Fig. S4), we purify the qubit state to 
1 ^) by measuring the qubit and applying a rotation if 
measured in |e). We succeed in preparing \g) with a proba¬ 
bility of 99%. Secondly, when performing qubit tomography 
we reset the qubit to \g) if it is measured to be in |e). Since 
we must wait Twait before feedback can be applied, the cavity 
state will acquire an additional phase Xqs^wd^it if the qubit 
is in |e). In this case, in addition to reseting the qubit, the 
FPGA applies an equivalent phase shift on the subsequent 
Wigner tomography pulse. This feedback implementation 
does not close the ‘locality’ loophole for a CHSH Bell test 
and therefore cannot be used to test local realism. 


2 Random and systematic errors 


Gaussian error statistics: We perform single-shot mea¬ 
surements that are discriminated into binary results and 
report measured observables taking the mean of N ex¬ 
perimental outcomes. All mean values are reported with 
N > 4000 measurements and the highest measurement fi¬ 
delity for any joint observable is = (V + l)/2 = 0.93 such 
that min[A^J^, A^(l — J^)] > 300 ^ 1. From this, we can 
approximate the mean value of all measured observables to 
follow a Gaussian distribution with standard deviation 

(Fig. S 6 ). 


Detector cross-talk: The sequential detection protocol in 
this experiment uses the same detector to perform first a 
qubit measurement followed by a cavity measurement. To 
minimize unwanted systematic errors due to detector cross¬ 
talk between measurements, we perform each experiment 
under four detector setting permutations. Two settings are 
used for the qubit measurement: a pre-rotation which maps 
a qubit eigenstate |±) to detector values and another 

which maps |±) to • Two settings are used for the cav¬ 
ity measurement: a Ramsey experiment which maps a cavity 
eigenstate |±) to detector values TMI and another which 
maps |±) to tTTI (Tab. S3). Each detector setting is per¬ 
formed an equal number of times and results are combined 
to remove unwanted correlations between detector readings 
and measured quantum observables. See Sec. 5 for an anal¬ 
ysis on the effects of these detector settings on a Bell test. 

The dominant form of cross-talk for this experiment is due 


to qubit state decoherence between measurements. To real¬ 
ize the cavity state measurement, the qubit must be initial¬ 
ized in 1^), which we perform using active feedback. Qubit 
decay can occur during this reset process causing an incor¬ 
rect initialisation for cavity state detection. We can model 
this error by observing the possible trajectories of each mea¬ 
surement outcome (Fig. S5). This modifies the average mea¬ 
surement of the observable AB where A, B are qubit and 
cavity operators that can be decomposed into qubit projec¬ 
tors AB = (A+ — A_)B, where A+ -f A_ = /. Due to 
qubit decay, the measured value {A_^B) will be modified to 
(1 — 2pc){A-^B) where Pc is the probability of qubit decay 
in the time between the first measurement and the feedback 
rotation. This relation changes the measurement into: 

^AB) ^(1 - 2pe) (A+5) - {A_B) (S 2 ) 

=(1 - Pc) {A^B - A_B) - Pc {A^B + A_B) 
=(1-Pe) {AB)-p, {B) 

For measuring B = Xc^Yc^Zc of the Bell-cat state 
we expect {B) = 0, which gives merely a reduction in the 
visibility of the observable (AB) by a factor(l —pc) without 
systematic offsets. We estimate in this experiment that Pc = 

_ ''"wait 

1 — e ^0.06. With this justification we can predict 

the additional loss in visibility V mentioned in the previous 
section which gives a visibility Vpred = ~ Pc) Y = 82%. 

The experimentally obtained visibility V is 85%; we believe 
the discrepancy between predicted and measured values is 
due to an overestimate in the time the qubit is susceptible 
to energy decay during measurement. 

Tomography rotation errors: We observe systematic ef¬ 
fects attributed to an amplitude error using the K!g opera¬ 
tion for pre-rotations used in qubit state tomography. A ar¬ 
bitrary pre-rotation = g^ cos 0 +^ 0 ; sm 0 ) transforms 

a qubit measurement along the Z axis: 

Z ^ cos 0Zsin 0 cos (l)Xsin 0 sin (/)¥ (S3) 

If ^ 7 ^ 0 or TT, a systematic offset can occur. This is ob¬ 
served in Fig. S7 where the contrast of {ZZc) is reduced 
and a residual offset in (ZXc) and (ZYc) respectively are 
produced. From Eq. S3, we predict the following relation¬ 
ship: 

(ZZc) tan e = / {ZXX + {ZY^ (S4) 

Erom measurements, we can approximate the fractional am¬ 
plitude rotation error 50 = 4.2%. This under-rotation is due 
to a photon-dependence of the calibrated pulse amplitude. 
Mitigating this error susceptibility is being explored with 
composite pulses in future experiments. 

3 state tomography 

Calculating observables: We can represent the density 
matrix in the excitation number basis: 

1 N 

p=Yl Y PiP 1 ^) ® 

i,j=0 n,m=0 
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where pfp are elements of qubit/cavity density matrix and 
|i,j) is the qubit state in the excitation basis and |n,m) 
is the cavity state in the excitation (photon number) basis. 
From a density operator, one can calculate an observable of 
the combined system by determining the product of observ¬ 
ables from each individual system: 

{AB) = Tr [ABp] (S 6 ) 

where ^4, B are operators for the qubit and cavity respec¬ 
tively. 

In the joint Wigner function, the qubit basis is the Pauli 
set ai = {/, }. For the cavity mode, we choose the 

displaced photon parity operator = D^PDl^ that corre¬ 
sponds to a single point in the cavity state Wigner function. 
For a truncated Hilbert space (in this experiment A^max = 12 
) and a displacement grid of (amax,min = ±3.4 with step size 
Aa = 0.085, this measured Wigner function represents an 
over-complete set of measurements for the cavity mode. The 
joint Wigner function Wi{a) = ‘^{(JiPa) is constructed di¬ 
rectly from experimental measurements. 

A qubit operator A can be written in the Pauli basis A = 
AiCTi where Ai = Ti[A(Ti] and a bounded cavity observ¬ 
able (see [38] for details) can be represented in continuous- 
variable basis B = ^ f B(a)Pad^a where B(a) = Pt[BPcx]. 
Finally, the composite qubit-cavity density matrix can be 
written as: 

P = ttY, [ Wi{a)aiPad‘^a (S7) 

i 

Note that for separable states p = Pq 0 Pc, this relation 
can be split up into their respective discrete and continuous 
components: 

Tr[/>gcri]cri (g) 27r / ±r[/)cPa]-Pad^a (S8) 

For any state p, we can write the mean value of an observable 
for the combined system with the following relation: 

{AB) = Tr [ABp] 

= Tr / AiB{a)Wj{a')aiCFjPaPa'd?ctd^ct' 

hj 

(S9) 

Using the following operator rules Tilaiaj] = 6ij and 
TilPaPa'] = S‘^{a — a') we can simplify Eq. S9: 

{AB) = E / AiB{a)Wi{a)d?a (SIO) 

i 

The overlap integral used in this calculation is similar to 
descriptions of the standard Wigner function [38, 9]. Shown 
in Fig. S7 is a comparison between observables calculated 
by Eq. S9 and those determined from a density matrix re¬ 
construction. 

Detector efficiency: Under experimental conditions, the 
measured joint Wigner function is determined with point- 
by-point measurements of the joint observable {oiPo). De¬ 
tector inefficiency results in a reduced visibility V G [ 0 , 1 ] 


and in turn a reduced contrast of the measured joint Wigner 
functions IU™®^®((a) = We can determine V by 

tracing over both the qubit and cavity states and comparing 
this to its ideal value f {a)d^a = 1: 

v = J Wf^^{a)d‘^a (Sll) 

where / is the qubit state identity operator. We observe 
V = 85% and attribute this primarily to readout infidelity 
and qubit decay between the sequential measurements (See 
Sec. 2 ). 

Density matrix reconstruction: In Eig. 2 of the main 
text, we show the reconstructed density matrix of a target 
Bell-cat state |'0 b)- We perform this reconstruction with 
a-priori assumptions that the cavity state is truncated to 
twelve occupied photon number states A^max = 12 , the re¬ 
sulting noise of each averaged measurement is Gaussian dis¬ 
tributed, and the reconstructed density matrix is positive 
semidefinite with trace equal to one. 

Under these constraints, we perform a least squares re¬ 
gression using a Maximum likelihood estimation [39]. To 
analyze this regression, we perform residual boostrapping 
on the reconstructed data set giving bounds on the error 
statistics of the inferred state. 

4 Encoded subspace 

Orthogonality of logical states: In the main text, we de¬ 
scribe encoded qubit states of the cavity where logical states 
|0l) , |1l) correspond to coherent states \f3) , \ —f3). This ap¬ 
proximation only holds for coherent states |±/d) that are 
quasi-orthogonal \{—P\P)\‘^ 1. To be more precise, we 

can calculate the maximum Von-Neumann entropy of the 
encoded space to determine its capacity to store informa¬ 
tion: 

^ [Pmax log 2 Pmax] 

= log2 Vi (S12) 

i 

where Pmax = \{\P) {P\ ± \~P) {~P\) is the density matrix 
for a complete mixture of the logical subspace and p is its 
set of eigenvalues. Rewriting Pmax in the even/odd cat state 
basis: 

Pma. = i(l + e-2|^l') |E) (E| + i(l - e-2|^l^) |0) (0| (S13) 

where |E), |0) = -^=i==(|/3) ± |-/3)). Recall that 

{E\0) = 0 for all coherent state amplitudes f3. This gives 
the following entropy relation: 

5 = -l(l+e-2|/5l^)log2 (i(l + e-21'^l')) (S14) 

_i(l_e-2|/5|=)lo&(i(l-e-2|/^l^)) 

Shown in Eig. S 8 is the capacity to store information using 
this encoding scheme. Entropy varies from zero bits to a 
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value asymptotically approaching a single bit with increas¬ 
ing coherent state amplitudes p. The orthogonality between 
logical states \{f3\ — /3)\‘^ is directly related to this informa¬ 
tion capacity and serves as a proxy for validating the qubit 
approximation of the produced cavity state. 

4.1 Encoded state observables 

The coherent state basis chosen in this report to represent 
the encoded qubit has Pauli operators: 

X, = \-i3){P\ + \f3){-f3\ (S15) 

Yc = j\-/3) (/3|-i|/3) (-/3| 

Zc = \/3) (/3I-I-/3) (-/3| 

Ic = \l3) {l3\ + \-l3) {-I3\ 

Here, we will show that the operators expressed in Eq. 4 
of the main text, approximate these encoded Pauli opera¬ 
tors. Assuming {f^\—/3) 1, we have the following photon- 

number parity P relations: 

(/3|Po|/3) = (/3|-/3)«l (S16) 

(/3|Po|-/3) = (^|^) = l 
(^1 1/3) = {I3- a|a - ^) « 1 

(/3|P«|-/3)=e2(“^*-“*« {a\-a) 


where Pa = D^PD^a for some displacement amplitude a. 
Now taking the projector M = \f3) {f3\-\-\—f3) (—/d|, we derive 
the encoded state’s Pauli Operators from the cavity state 
observables reported in Eq. 4 of the main text: 

MPoMt « 1-^) (^1 + 1/3) (-/3| (S17) 

MP^Mt « 1^) (/3| 

MP.^sMt « I-/3) (-/3| 
MP^Mt«j|-^) (/3|-j|/3) (-/3| 
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Putting these relationships together, as in Eq. 4, builds the 
encoded state observables {Xc^Yc, Zc, Ic} and reveals that 
these observables can be efficiently measured using Wigner 
tomography. Ic and Zc require a comparison between two 
different observables. Eor true single-shot readout of these 
logical observable Zc, measuring a single value in the cavity 
state Husimi-Q distribution Q{/3) = ^ {/3\p\/3) can be em¬ 
ployed where Zc = 27rQ{/3) — 1 . This is being explored in 
future experiments. 

Encoded state Pauli set: We can represent the two qubit 
Bell state shown in Eig. 2 as a list of two-qubit correlations. 
The complete set constitutes the permutation of each of 
the single qubit Pauli set {/, X, T, Z}. We can determine 
the two-qubit Pauli set from the complete reconstructed 
qubit-cavity state and projecting onto the encoded basis 
of {Ic^Xc^Yc, Zc}. Eig. S7 shows the resulting two-qubit 
Pauli set for the transmon qubit and an encoded qubit in 
the cavity mode, a variant of the reduced density matrix 
representation shown in Eig. 2 of the main text. 


Encoded state preparation: We can diagnose errors that 
can occur during state preparation from the reconstructed 
Pauli set. The dominant nonideal effects we explore are 
qubit decay during preparation and single-qubit rotation er¬ 
ror. 

During state preparation, the product state \ 2 p) = 

^(1^) -T |e)) (8) \/3) is initialized. Under the dispersive 
interaction, the system evolves into the entangled state 
To describe the effects of Ti decay, 
we can look at the diagonal elements of the reduced den¬ 
sity matrix after the entangling evolution of the dispersive 
interaction: 

diag[p] = |{ \g,l3){g,/3\ + e“^ |e,-/3) (e,-/3| (S18) 

+ 'ZCk \g,ak) {g,ak \ } 

k 

where ai = represents coherent states when a jump 

occurred at time tk and 7 = Projecting onto the logical 
basis (here we will approximate \ {ak\P)\‘^ ^ 1) produces the 
resulting scaling on the joint Pauli measurements for the 
ideal state: 

(lie) (X i (1 + e-^) (ZI,) (X i (1 - e-^) 

(S19) 

(ZZe) (X pi + e-^) {IZ^) X pi - e-^) 

This gives us an approximate method to predict our abil¬ 
ity to prepare a state that is within the logical subspace 
given our experimental parameters, {He) = 0.99. Our mea¬ 
sured value taking into account detector inefficiencies pro¬ 
duces {lie) = 0.98. 

The preparation of this entangled system is also sensitive 

2L 

to the amplitude of the initial qubit rotation Y/2 = . 

The angle of rotation 0 will determine the prepared state 
as: 

I//)) = A/'{cos I \g, /3) + sin f |e, -/3)} (S20) 

Eor the states prepared nearly as I'^b), ^ ~ f (1 “ 

will result in the following modification of the joint Pauli 

measurements for an ideal state: 

{ZIe) cc 160 

(S 21 ) 

(/Ze) X ^59 

Erom the measurements in Eig. S7, we can determine that 
the relative error for the rotation angle in our preparation 
rotation to be SO ^ 2 . 8 %. 

5 Bell test analysis 

The main text reports CHSH Bell tests composed of two 
qubit observables A, B and two cavity observables Ac, Bc^ 
correlated such that: 

O = {AAc) + {ABc) - {BAc) + {BBc) (S 22 ) 

We perform two variants of this test on the state |'0 b)- 
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Test Model: In the first test we choose qubit cavity ideal case is the result of four summed joint Wigner values 
observables and qubit observables Z{0),X{0) where: represented as: 


Z{0) = Z cos ^ — X sin | X{0) = X cos | ^ sin | 

(S23) 


Oideai = 2 (cos4q!o/ 3 + sin4Q;o/3)e 


This angle 0 corresponds to a rotation of the qubit state 
before detection. In Fig. 4a, we plot O for each of the four 
permutations of the joint observables and find a maximum 
Bell violation for an angle 0 = giving observables: 


x-z 


A — X-\-Z . TD — _ 

“ x/2 ’ ^ ~ V2 

Ar — Zc\ Bp = Xr 


(S24) 


where is an optimal displacement for maximum viola¬ 
tion which can be calculated from Eq. S27 and in detail in 
Ref. [21]. Taking into account photon loss and measurement 
inefficiency produces the following relationship: 

Opred = (cos4(ao/d + sin4(ao/d)e“^l^°l 


As shown in Fig. 4c, we can model the effects of photon 
loss and measurement inefficiency on the maximum viola¬ 
tion. For the ideal case, an overlap of the coherent state 
superposition reduces contrast in {AZc) and {BZc) and will 
limit the maximum Bell signal: 

©ideal = \/2(2 - 

Measurement inefficiency will reduce the contrast of this 
maximum Bell signal which we expect to go as the visibility 
V: 

©vis = V2V(2-e-®l^l') 

Photon loss will also have an effect on the maximum Bell 
signal by reducing the measured contrast of all correlations 
for {AXc) and (BXc). This produces the an amplitude de¬ 
pendent maximum Bell Signal: 


This predicted Bell signal is shown in Fig. 4b of the main 
text using the measured joint-Wigner contrast V = 0.85 and 
an effective time teff = 1.24 /rs. 

Optimal measurements for encoded observables: 

Eq. 3 of the main text describes the ideal observables to 
efficiently determine an encoded qubit state observable us¬ 
ing a superposition state with \/3\ ^ 1. In fact, the optimal 
measurement for particular observables will be further mod¬ 
ified for smaller coherent displacements. 

For the second CHSH experiment, the optimal observable 
B±jao ^ A ^c) follows the relation: 

^ —^=tan4(ao/d (S27) 

P + Oq 


where 7 = ^ such that is the photon decay time constant 
and teff is the effective time to create and measure the Bell- 
cat state. Finally taking into account both visibility and 
photon loss produces the expected maximum Bell signal: 

©pred = V2V(1 - - e-2T'l/5l') 


where is the amplitude for a coherent displacement Djo^^ 
to perform the measurement Pjao given p. Further details 
are discussed in Ref. [ 21 ]. In the large [3 limit, the observable 
corresponds to the encoded qubit state observable ^ {Xc + 

Yc) and follows the relationship as related in Eq. 4 of 
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the main text. Shown in Fig. S9 is the predicted and chosen 
optimal values for a maximum CHSH Bell signal. 


This predicted Bell signal is shown in Fig. 4 of the main text 
using the measured joint-Wigner contrast V = 0.85 and time 
between cavity state creation and detection teff = 1.24 ps. 

Test #2 Model: In the second test, we choose qubit ob¬ 
servables X, T and cavity observables Xc{a),Yc{a) where: 

Xc{a) = DjaPoD]^ « Xc cos ^ + V sin ^ (S25) 

Yc{a) = DjccP^D]^ « Yc cos - V sin ^ 

Where the displacement amplitude a corresponds to an ap¬ 
proximate rotation of the encoded cavity state before detec¬ 
tion. In Fig. 4b, we plot O for each of the four permutations 
of the joint observables and find a maximum Bell violation 
for a displacement a = 0.15for/d = l which produces the 
approximate observables: 


A = X- 

A _ X,+Y,. 


B = Y 

Br = 


Y2 


(S26) 


Shown in Fig. 4c, we can also model the effects of photon 
loss and measurement inefficiency for the second test. The 


Two-qubit entanglement witnesses: Two qubit entan¬ 
glement can also be quantified by an entanglement wit¬ 
ness W = lie — XXc Y YYc — ZZc [40] for a Bell state 
I'lp) = -^{\gg) + The witness ‘confirms’ entanglement 

for all observations of (W) < 0. Shown in Fig. SIO, we report 
W (as well as its corresponding direct fidelity estimation P) 
as a function of coherent state amplitude [3 using the optimal 
displacements described in Fig. S9. As expected, entangle¬ 
ment is not detected for a yd = 0 coherent state (a product 
state ■^(\g) + |e)) ® | 0 )). 

Bell test for each detector setting: We analyze the 
systematic errors that can occur from a particular detector 
setting. Shown in Fig. 5 are the observables used to cal¬ 
culated a Bell violation using test 7^2 for each of the four 
detector settings Sec. 2 . Systematic errors are shown to be 
within stastical bounds of the experiment and each detector 
setting violates Bell’s inequality by at least three standard 
deviations, see Fig. 5. In the main text, we report mea¬ 
surements from the combined data set resulting in smaller 
statistical error and a stronger violation of Bell’s inequality. 
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6 Quantum measurement back- 
action 


The sequential measurement protocol allows us to observe 
the result of quantum measurement back-action of the qubit 
on the cavity state. The result of an ideal qubit measure¬ 
ment outcome Mm will give a projected qubit-cavity state: 


I / \ Mm I'ip) 

iV^m) = ! = 


(S28) 


Measuring along the {X, T, axes of the qubit gives three 
measurement sets: 


X : 


Y : 


Z : 


1 1 
1 1 


1 
3 

1 0 
0 0 


-3 


1 -1 

-1 1 


1 ( ^ ^ 
2\-j 1 


0 0 
0 1 


(S29) 


Bell-cat projections: We prepare the system in a Bell- 
cat state as in Eq.2, and measure along each of the three 
qubit axes. These three measurements results in six possible 
outcomes |'0rn) = |'0g)^|'0c) with the projected cavity states: 


iV^cav) X \ N (1/3) -h 1/3)) 
F: V (1/3)-j 1/3)) 
Z: 1/3) 


V(|^)-|/3)) 

V(|/3)+i|/3)) 

I-/3) 


(S30) 


See Fig. S13 for each projective measurement of the Bell-cat 
state |'0 b)- The method of using strong projective measure¬ 
ments to create cat states has been demonstrated in previous 
experiments [41]. 

Fock state projections: We prepare the system in a state 
such that the qubit state |e) is correlated with the pho¬ 
ton Fock state \m) of a coherent state \P) (in this example 
m = 3 photons and P = v^). This can be written as: 


IV’) = Cm |e, m)+'Y^Cn \9, n) (S31) 

n^m 


where Cm = (^|/^)- Shown in Fig. S14, when the qubit is 
measured along the Z axis we observe a change in photon 
statistics such that a +1 event projects the cavity onto the 
state IV^cav) = A/’dyd) — Cm 1^)) and a —1 event projects 
onto the Fock state |'0cav) = |^)- 
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Figure SI: Photograph of device and amplifier: (a) One half of the 3D circuit QED device shows both the readout 
and storage cavities. Strongly coupled to each cavity is a single vertical transmon. (b) High-fidelity measurements are 
achieved with near quantum-limited amplification provided by a Josephson bifurcation amplifier. Shown is the chip and 
sample holder for this device. 



Figure S2: Experiment Schematic 
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Figure S3: Feedback inset: The feedback setup uses two input-output (I/O) boards for qubit and storage resonator 
control and one arbitrary waveform generator (AWG) for readout resonator control. All have a dedicated microwave 
generator and mixer for amplitude and phase modulation. Each I/O board has five main components: 1) a digital- 
to-analog converter (DAC) for pulse generation; 2) digital outputs serving as marker channels; 3) an analog-to-digital 
converter (ADC) that samples input signals; 4) an FPGA that demodulates the signals from the ADC and based on 
predefined thresholds determines the measured qubit state, \g) or |e) to generate pulses; and 5) a PCIe connection that 
transfers FPGA data to a computer (PC) for analysis. In this setup the top I/O board serves as the master, which accepts 
the readout signal, returns qubit state information, and using digital output signals, triggers the AWG and the second 
I/O card given a particular qubit measurement result. 


Term Measured 
(Prediction) 
ujq/27r 5.7651 GHz 

ujs/27r 7.2164 GHz 

ujr/27r 8.1740 GHz 


Kq/27r 240 MHz 



1.5 kHz 


Qubit 

Storage 

Readout 


(2 kHz) 

Ti 

10/is 

- 

- 

X9s/27r 

1.43 MHz 

T2 

10/is 

- 

- 

Xgr/27r 

1 MHz 

T:av 

- 

55/is 

30ns 

Xrs/27r 

(1.7 kHz) 

ground state (%) 

90% 

> 98% 

> 99.8% 


Table SI: Hamiltonian parameters Table S2: Coherence and thermal properties 
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Figure S4: Pull quantum circuit diagram Each experiment is split into four components, (a) First, the system is 
initialized. The qubit state is measured and a qubit pulse RS = {R^ or 1} is applied to reset the qubit to \g). (b) Second, 

the entangled state is created with a cavity displacement and a qubit rotation Yl2 = R^J^ followed by a ^ waiting 
time to produce the entangled state I^/^b) = (c) Following preparation, a qubit state detection is 

performed with a pre-rotation Ri (see Table S3), a measurement, and a qubit reset RS. Finally, we perform a cavity state 
measurement using Ramsey interferometry where Rj = ±T/2 combined with an initial pre-displacement Dc^. This maps 
Pa to the qubit state which is read out with a subsequent qubit measurement. Correlations are reported as the product 
of detection events between measurements in (c) and (d). 


p 



Figure S5: Measurement trajectories given qubit decoherence. We can model the behaviour of qubit decoherence 
in a single measurement trajectory. Qubit decay (which occurs with a probability p) can lead to an improper initialisation 
of the second detection and in turn produces an incorrect measurement result. This form of detector cross-talk can lead to 
a reduction in visibility and potential systematic offset of the measured qubit-cavity observable: {AB) (1 — Pc) {AB) — 
Pc {B) (see section 2). 
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Residuals from reconstruction {SW^{a)) 


Qubit Cavity 



Figure S6: Histogram of reconstruction residuals. 

Plotted are the residuals corresponding to the density ma¬ 
trix reconstruction of the Bell cat state shown in Fig. 2 of 
the main text. This Histogram shows the distribution of the 
25, 000 residuals from the joint Wigner function which gives 
a Gaussian distribution (mean value /i = 7.0 x 10“^, stan¬ 
dard deviation a = 0.015), which agree with our expectation 
for statistical error (Test = ~ 0.015. 


Table S3: Table of measurement operators. 

As shown in Fig. S4, pre-rotations before qubit 
and cavity state measurements determine the mea¬ 
sured observable. Shown are the different pre¬ 
rotations used and the corresponding measure¬ 
ment operator. 
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Figure S7: Reconstructed Pauli sets. The set of sixteen joint Pauli operators span the two-qubit Hilbert space of the 
qubit/encoded-qubit state. Shown is the Pauli set for the entangled target state derived in two ways. (Red) is the 
reconstructed Pauli set using a density matrix reconstruction of the full quantum state with no normalization constraint, 
then projecting onto the encoded subspace. (Blue) shows the values discerned from an overlap integral of the measured 
joint-Wigner functions (Eq. S9). These measurements agree with each other within statistical errors. 
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Figure S8: Informational entropy. The capacity to store 
information into a cat state is determined by the orthog¬ 
onality of its logical states |/3), |—/d). Shown is a compar¬ 
ison between the coherent state overlap (dashed line) and 
the maximum Von Neumann entropy Eq. S14 (solid line) 
for this logical encoding. Notice that entropy rapidly ap¬ 
proaches one bit for > 1, ensuring that information can 
be reliably encoded into the coherent states with manage¬ 
able separations. 


Figure S9: Optimal displacement for Bell violation. 

For performing Bell test 7 ^ 2 , the optimal observables to 
measure maximum correlations depend on the size of the 
Bell-cat state Eq. S27. The dashed black line shows numer¬ 
ically calculated optimal displacement points as a function 
of coherent state amplitude /3. Shown in circles are the ex¬ 
perimentally determined optimal displacement values used 
to measure a maximum Bell violation. Differences between 
chosen and ideal values are a result of the discretization of 
our measurement settings. The dashed green line is the ap¬ 
proximate trend |(aopt| = |^| foi* large cat states, which 
diverge at small [3. 
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Figure SIO: Entanglement witnesses with cat states. An entanglement witness and direct fidelity estimation (DFE) 
are determined by measuring four qubit-cavity correlations, (a) The entanglement witness W = II — ZZ — XX YY 
shows entanglement for all negative values (grey shading), (b) DFE to a target Bell state T = II Y XX — YY + ZZ is 
also shown where entanglement can be confirmed for values above T > 0.5. Notice that these two witnesses have a much 
looser bound for entanglement than the CHSH Bell test. 
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Figure Sll: Observables from each detector setting. 

To ensure that a particular detector setting is not producing 
systematic errors we have not taken into account. We report 
a Bell test for each detector setting used to observe our 
maximum violation in test ^2. The expectation value of 
each observable used in that Bell test is shown for the four 
detector settings used Sec. 2. Significant deviations due to 
unexpected systematic errors are not observed. 


Figure S12: Bell test for each detector setting. A Bell 
test is analyzed for each detector setting to determine the 
effects of possible systematic errors. Each of these subtests 
violate Bell’s inequality by more than three standard devi¬ 
ations of their statistical error. 
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Figure S13: Qubit measurement back-action of a Bell-cat state. The resulting projections of the state I^Pb) = 
/^) + ~P)) fo a particular qubit measurement outcome. Note that measuring along the X and Y axes results 

in a projected cat state each with different superposition phases. Combining these measurements with the probability to 
obtain each result describes the entire system and is used to create the joint Wigner function representation in Fig. 2 of 
the main text. 
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Figure S14: Qubit measurement back-action for an entangled Fock state, (a) A measured Wigner function of a 
coherent state \/3) where P = \/3 results in a Poissonian photon distribution. Performing a photon-selective qubit rotation 
on the level where m = 3 results in an entangled state j^/;) = Cm |e, m) + ^n^m \9^ p) where Cn is the coefficient 
of the photon number state Cn = (^|/d). (b) The measured Wigner function of the cavity state after the qubit has 
been measured in the —Z state results in a 3-photon Fock state, (c) Instead, when a result is obtained the measured 
cavity state Wigner function is a Fock-state subtracted coherent state |^)- 
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